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SUMMARY 

Control laws were developed to investigate methods of alleviating the response of a 
STOL transport to gusty air. The aircraft was simulated in a low-speed approach con- 
figuration. The elevator, spoilers and trailing edge flap element were used as control 
devices. Two accelerometers, one mounted vertically and one mounted longitudinally 
with respect to the airframe, pitch angle and pitch rate gyros, and an airspeed sensor 
were used to generate feedback control commands. A conjugate gradient optimization 
procedure was used to select the feedback gain matrix which minimized the average of a 
quadratic functional involving normal acceleration (at both the C.G. and tail of the air- 
craft), longitudinal acceleration, pitch angle and rate, inertial speed, and glide path 
angle variations. The optimization was performed for a selected gust field consisting of 
vertical and horizontal gusts in the form of series of constant amplitude pulses with 
various durations. 

Feedback matrices for various combinations of elevator, flap, and spoilers were 
obtained. The root-mean-square (rms) magnitude of several variables appearing in the 
performance index are presented for each control configuration. The closed-loop pole 
locations and magnitude versus frequency plots of the aircraft transfer functions relative 
to horizontal and vertical gust inputs were obtained and compared. 

Feedback matrices were also obtained for a number of sensor configurations uti- 
lizing pitch angle and rate and one or more of the following variables: normal accelera- 
tion, axial acceleration, relative airspeed at the center of gravity or at the nose of the 
airframe, and angle of attack at the nose of the airframe. The closed loop system 
characteristics for these sensor configurations are compared with the case utilizing only 
pitch angle and rate feedback. 

The study was performed on the EAI 690 hybrid computer facility at the Langley 
Research Center. The linearized equations for the longitudinal dynamics of the aircraft 
and die control system were simulated on an analog computer, and the optimization 
algorithm was programmed on a digital computer. The modal and frequency response 
computations w^re performed on a large general purpose digital computer. 
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INTRODUCTION 

The reduced wing loading required for STOL capability increases the aircraft’s 
sensitivity to atmospheric gust disturbances. This, in turn, increases the motions and 
accelerations imposed upon the passenger compartment and adds to the level of pilot 
work load in gusty weather. In addition, the low speed of operation results in a reduced 
aerodynamic control surface effectiveness. In view of these factors, a study was made 
of the longitudinal control of a STOL aircraft. The object of this study was to determine 
the benefits to be realized by utilizing aerodynamic control surfaces on the wing and tail 
to improve the stability and ride quality. This study considered the reduction of air- 
craft motions relative to a steady flight condition and did not address the problem of 
flight control by the pilot. 

The geometric characteristics of the high-wing STOL aircraft examined in this 
study are presented in figure 1. The aircraft has triple-slotted, externally -blown jet 
flaps and a large T-tail. A detail of the wing showing the configuration of the engines, 
spoiler, and flap assembly is presented in figure 2. The three control devices consid- 
ered in this study were thr» spoiler, the trailing -edge element of the flaps (hereafter 
referred to as the flaps), and the elevator. The spoiler provided, primarily, lift modula- 
tion. The flap element provided substantial modulation of both lift and drag. The eleva- 
tor acted as a conventional rotation controller. The engines were fixed at the trim thrust 
level. 

Aircraft simulation and control system configuration studies were performed on 
the EAI 690 hybrid computing facility at Langley Research Center. The linearized longi- 
tudinal dynamics of the transport, equipped with fast acting elevator, spoilers, and flaps, 
were programmed on the analog portion of the hybrid computer (EAI 680). The aircraft 
was configured for landing approach. First order lags were used to represent actuator 
dynamics. A feedback control matrix operated on a number of sensor signals to produce 
control command signals. Sensor variables, selected flight variables, and control com- 
mand variables were tabulated at regu’ar intervals by the digital portion of the hybrid 
computer (EAI 640). Two series of rectangular test wave forms for horizontal and verti- 
cal gust inputs were generated according to stored tables in the digital computer. A 
quadratic performance index and its gradient were evaluated, off-line, by the digital 
computer. A conjugate gradient optimization procedure (ref. 1) was used to generate 
feedback gain matrices minimizing the performance index. 

Feedback gain matrices were computed for the cases of elevator, spoilers, and 
flaps each operating alone, in combinations of two, and all three together. The magni- 
tude of each of the variables contributing to the performance index for each combination 
of control devices is presented. Closed loop pole locations and magnitude versus fre- 
quency plots of the aircraft transfer function relative to sinusoidal horizontal and vertical 
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gusts are examined for each feedback gain matrix obtained above. These items are i 

compared to determine the relative effectiveness of each control configuration in modi- j 

fying the gust response of the aircraft. 

Feedback gain matrices were then computed for elevator, spoiler, and flap controls 

V 

and several sensor configurations. The case utilizing pitch angle and rate variables as 
feedback quantities was used as a baseline for comparison. Signals from normal and 
‘ axial accelerometers, airspeed sensors at the nose of the airframe and at the center of 

\ gravity, and an angle of attack sensor at the nose of the airframe were added individually 

>' to the basic sensor configuration. In addition, both the airspeed sensor at the nose and 1 

the normal accelerometer and then the airspeed sensor at the nose and the angle of attack i 

sensor were added to the basic configuration. These cases were compared to determine j 

the effects of each sensor addition upon the stability of the system and the system 
response to gust inputs. i 

. s 
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SYMBOLS | 


a 0 normal acceleration at rear passenger station, g units 


a n normal acceleration at the C.G., g units 


axial deceleration, g units 


c mean aerodynamic chord, meters (feet) 


Cr. drag coefficient, — 

D qS 

C T lift coefficient, — 

L qS 

C pitching moment coefficient, 

m qcS 

aC 

C . * — example of aerodynamic derivative 
m e ae 

C engine thrust coefficient 

r* 


C.G. center of gravity 


D 
g 

~ v 


drag force, newtons (pounds) 

o o 

gravitational acceleration, m/s (ft/sec ) 



gradient vector for ith approximation of the solution set of gains 
feedback gain matrix 

element of the gain matrix identified by two indices 

ith approximation of the solution set of gains 

vertical distance from C.G. to the horizontal tail divided by c 

incidence of tail, degrees 

performance index 

index denoting sample time of variable 
radius of gyration divided by c 

lift force perpendicular to velocity vector, newtons (pounds) 
distance from the nose of the aircraft to the C.G. divided by c 
distance from C.G. to rear passenger station divided by c 
longitudinal distance from C.G. to aerodynamic center of tail divided by c 
pitching moment, N-m (lb-ft) 
mass, kg (slugs) 

number of samples used in evaluation of performance index 
dynamic pressure, N/m (lb/ft ) 
symmetric weighting matrix on control commands 
root mean square value 

symmetric weighting matrix on measured variables 


i 


ith search direction in optimization algorithm 

2 2 

area of the tail, meters (ft ) 

2 2 

area' the wing, meters (ft ) 
time 

time constant of elevator actuator, seconds 

time constant of flap actuator, seconds 

time constant of spoiler actuator, seconds 

nondimensional speed variation, AV/V 0 

nondimensional airspeed measured at C.G., V^/V Q 

nondimensional airspeed measured at nose of aircraft, Vp/V Q 

stream velocity variation at tail due to gusts, m/sec (ft/sec) 

nondimensional horizontal gust component, V V o 

nondimensional vertical gust component, Vy/V Q 

vector of control commands 

inertial velocity, m/sec (ft/sec) 

airspeed measured at C.G. 

airspeed measured at nose of aircraft 

horizontal gust component 

vertical gust component 

inertial velocity variation, m/sec (ft/sec) 


vector of performance variables 


vector of desired performance variables 
vector of sensor signals 
angle of attack, radians 

portion of angle of attack due to aircraft motion, radians 
angle of attack variation, radians 

portion of angle of attack variation due to aircraft motion, radians 

sensed angle of attack at nose of airframe, radians 

flight path angle, radians 

flight path angle variation, radians 

elevator variation, radians 

flap variation, radians 

spoiler variation, radians 

downwash angle, radians 

weighting on gain matrix elements in performance index 

pitch angle, radians 

pitch angle variation, radians 

air density, kg/m^ (slugs/ft^) 

one -dimensional search parameter 


defined as i t c/V Q 



Subscripts: 


c control command signals 

i,j,£,m indices used in identifying gain matrix elements 

o nominal flight variables 

t tail parameters 

w wing (tail-off) parameters 

Superscripts: 

d desired flight variables 

o nominal parameters and variables in optimization algorithm 

T transpose of a vector or matrix 

' (prime) vector evaluated for a perturbed set of gains 

A dot over a variable denotes the first derivative with respect to time. 

Two dots over a variable denotes the second derivative with respect to time. 

SIMULATION EQUATIONS 
Aircraft Model 

The mathematical model of the aircraft used in this study is developed in appen- 
dix A. The longitudinal dynamics are linearized about a steady flight condition and are 
referenced to the wind axis coordinate system (see fig. 3). The flow variations in the 
vicinity of the tail (downwash and stream velocity) are considered to be delayed by the 
time required for a change occurring at the wing takes to have an effect at the tail. These 
delays are modeled at first order lags. The equations of the aircraft model follow. 
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The actuator dynamics are modeled as first order lags. The dynamic equations 
for the control surface deflections are 


« ■ («c - 5t )/‘f 

Ss = (Ss c - 5s)/t s 

«4 . (ee c - «.)/», 


where 6f , 6s , and 6e are the flap, spoiler, and elevator command signals, 

V v C 

respectively. These signals are to be generated according to a feedback control law. 
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Sensor Signals 

Pitch angle and pitch rate gyros are assumed to provide measurements of o and 
0. An accelerometer is mounted vertically with respect to the fuselage at the aircraft's 
center of gravity. The signal produced is a mixture of acceleration components in the 
lift and drag directions: 


a 

« g 


^0 - dtj^cos a 0 - u sin ajj 


where a R is measured in g units. 

A second accelerometer is mounted longitudinally with respect to the fuselage at 
the aircraft center of gravity. The signal produced by this sensor is: 


*X = g [(® • "l) 8in a o + “ cos -.] 


where a^ is measured in g units. 

An airspeed sensor is placed at the center of gravity. The sensed signal is the 
magnitude of the velocity of the air relative to the aircraft: 

U A = u + U H cos ‘ U V sin y o 

Relative airspeed and angle of attack sensors were placed on the nose of the air- 
craft. These sensors provide measurements of gust inputs before they effect the forces 
on the aircraft. The equations describing the outputs of these sensors are 


Up = u + u H cos - u v sin y n - 


n d 


o “V 


sin a. 


6/ n c 


a F = + U H 8in y o + U V cos y o " ~T~ cos “o 


All sensors are assumed to be ideal with infinite bandpass and no lags. 
The inertial flight path angle is 


r= r o+ 0 


°i 


The normal acceleration of a point in the fuselage located l wing chords behind 

r 


the center of gravity is 


a B = a n ' 


l p cQ 
g 


i 



Control Law 


A linear control law is used to generate the actuator signals for each control con- 
figuration investigated. The general form of the control law is 


u = Gz 


where 

u vector of actuator command signals 
G control gain matrix 
z vector of sensor variables 

The elements of u and z and the dimensions of G are selected as appropriate for 
each control configuration. 


Performance Index 


The quadratic performance index used in this study has the general form 


where 



vector of all measured variables at sample instant k 

y£ vector of desired values for the measured variables at sample instant k 
u k vector of control commands at sample instant k 
G^j ijtli element of the feedback gain matrix, G 

S,Q symmetric weighting matrices on the state errors and control commands, 
respectively 

tj weighting on the gain elements 


and 


l 


!» 

e 

e 

U A 

y 

u 

a B 


u * 
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The performance index is composed of terms weighing the difference between the 
desired values of the measured variables and their actual values, anc the values of the 
control command signals with a penalty placed on the magnitudes of the gain elements. 
The penalty restricts the magnitudes of the elements of the gain matrix to acceptable 
values. In the absence of this term several gain elements have invariably taken on 
excessive values with little change in the control commands. 


Given a small perturbation in the value of the gain element G^ m from its normal 
value G° m , the measured variables and the control commands are assumed to attain 
new values according to 


Zk * & * s X k ■ z£ - 


“k = Hk ♦ s Hk * “k 
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aG im 

’ du " 


6G. 


im 


'im 


6G. 


im 


where the superscript o denotes variables evaluated for a nominal gain matrix G°. 


Given this assumption, the first order variation of the performance index due to a 
perturbation of the gain element, G^ m , can be expressed as 

... 

n k -1 


6«T 


The term 


tel 


IZk - z£) T s 


% 


dG lm 


(u°) T Q 


dU 


dG i 


m 


, + t)G. 


’fm 


6G. 


im 


aj c 


?G 


im 


6G, 


'im 


is the variation of the performance index J with respect to the 


gain element G^ m about G = G . 
The matrix 



"aj° 

8J° 


*7l 

aj e 


K * 
*0 

Wi 

• 



• 

9 



is defined as the gradient of J with respect to the matrix G evaluated at G = G°. 
This is analogous to the definition of the gradient with respect to a vector. 
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The value of an element of g Q 
element of G by a small amount and 


is approximated by changing the corresponding 
computing: 


aj 

aG fm 




- s) ♦ (sM* - 4 


aG, 


7m 


+ 


where yj^ and u^ are the vectors y^ and u^ evaluated for the gain matrix G° 
with the element G° m replaced by G° m + CG^. 

Notice that neither the evaluation of the performance index nor its gradient requires 
explicit knowledge of the plant parameters. This technique could be used in simulations 
utilizing hardware elements without requiring explicit mathematical models of the 
elements. 


COMPUTATIONAL TECHNIQUE 

A conjugate gradfent method of minimizing a performance index with respec to the 
elements of a feedback matrix was used to generate feedback matrices for several con- 
trol and sensor configurations. This method is attractive since it exhibits second order 
convergence with the amount of computation comparable to first order methods. The 
method used in this study was patterned after the method of reference 1. It is a first 
order algorithm similar to s'eepest descent in that successive search directions are 
based on gradient of the performance index at each step. 

The simple steepest descent methc l uses the current gradient as the search 
direction: 


s i = ^i 

where g . is the gradient of the performance index with respect to the gain elements for 
the ith gain matrix and Sj is the ith search direction. 

The conjugate gradient method of reference 1 improves on the steepest descent 
method by utilizing information contained in the previous search direction in forming a 
new search direction. The algorithm used in the study computed the search directions 
according to: 

S 1 = -&1 

( g i’ g i) 

>l ” gl + («i.r > i -i ) ai - 1 
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where (g^,gj denotes the inner product of with itself and (gj_i.Sj.jj denotes the 
inner product of g._ j with s._j. A single parameter minimization is performed along 
the search direction at each iteration step. The solution of this search is the nominal set 
of gain elements for the following step. 

The computational logic flow chart for implementing this conjugate gradient algo- 
rithm is given in figure 4. This chart includes some reasonableness tests and restart 
loops that are not discussed here. 

AIRCRAFT PARAMETERS 

The aircraft parameters were based on references 2 and 3. The configuration has 
full span, triple slotted flaps with the trailing flap element (hereafter referred to as the 
flap unless otherwise specified) used as a control device. Spoilers are mounted ahead of 
the flaps to provide direct lift control. Aerodynamic coefficients were obtained for a 
nominal flap setting of 25°/10°/50° for the three flap elements, engine thrust coefficient 
C ljL = 1.87, spoiler prop-up (up- rig) of 15°, and angle of attack of 10°. This condition 
results in a 3.88° glide slope at 35.4 m/sec for a 25 000 kg aircraft. Estimates of down- 
wash angle characteristics and of tail plane and elevator lift and drag were empirically 
derived from references 2 and 3. A complete list of parameters for the study case is 
presented in appendix B. 


COMPUTER OPERATIONS 

The simulation of the aircraft dynamics, actuators, and sensors was performed on 
an EAI 680 analog computer. This was interfaced with an EAI 640 digital computer 
through digital -to -analog and analog -to -digital channels and binary control lines. Fig- 
ure 5 illustrates the configuration of the hybrid computer system for automatic simula- 
tion and control law optimization. The aircraft dynamics, sensors, and actuators are 
simulated in the conventional analog fashion. Gust inputs are generated by the digital 
computer. Control loops from a number of sensors to the control actuators are imple- 
mented by digital -analog multipliers. These devices are used as analog attenuators 
whose gains are manipulated by the digital computer. Sensor outputs and other vari- 
ables used in the performance index are routed to the digital computer through analog - 
to-digital converters. 

A setup and static check program was written in the HYTRAN Operations Inter- 
preter Language (HOI). This digital computer program processes the aerodynamic 
coefficients, derivatives, and kinematic data to determine the potentiometer settings for 
the analog simulation. The potentiometers are set to these values and series of static 
checks of the analog setup are made. Test signals are applied to the initial condition 
terminals of each integrator and to each input to the simulation. The output of each 
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potentiometer and amplifier and the input to each integrator is checked by the program. 
This procedure identifies each component whose output is out of tolerance or which is 
wired incorrectly. 

A simulation control and optimization program was written in FORTRAN. Upon 
satisfaction of the static checks, this program is compiled, loaded, and executed. This 
program provides outer loop control of the analog simulation, generates the test gust 
signals, evaluates the performance index and its gradient, and varies the feedback gains 
according to the conjugate gradient optimization procedure. 

The simulation of the aircraft was performed for a flight interval cf 20 seconds 
duration. This is subdivided into 100 subintervals of 0.2 second each. The test gust 
signals were set before each subinterval. Then, the sensed variables were transferred 
to the digital computer through analog-to-digital converters. The digital computer com- 
mands the analog to resume simulation. Upon completion of the 0.2 second subinterval, 
the analog goes into the hold mode and sends a signal to the digital computer. 

After commanding the analog to simulate, the digital computer stores the sensed 
variables and computes their contribution to the performance index. Then the digital 
computer waits for the interval completion signal from the analog before obtaining the 
next set of sensed variables. After 100 such intervals, the stored performance informa- 
tion is used in the optimization routine. The optimization routine continues until manu- 
ally terminated. 


RESULTS AND DISCUSSION 


Effects of Control Configuration 

A series of optimization runs were performed for various combinations of flap, 
spoiler, and elevator to determine the effects of each control combination on the perfor- 
mance of the optimized systems. The performance index used for this series was 


■™ k £*{ 200 a "( k ) 


+ 200 e 2 (k) + 100 0 2 {k) + 100 a 2 


x (k) 


+ 100 y (k) + 200 




100 a i + 100 6f? + 20 fie? 


B 


(k) 


'00 


'(k) 


+ 20 6s 2 > + 0.002 Z (Gjj) 2 


: (k) 


i* j 


where the (k) refers to the values of the variables at the 'k+l) st sample time. This index 
is comprised of terms relating to ride qualities, a n > 0, a^, and ag; flight parameters, 
0, y, and V; control command signals, fie c , 6f c , and 5s c ; and the magnitudes of the 
elements of the gain matrix, Gy. 
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The simulation of the aircraft was initialized in steady, trimmed flight. Two pre- 
determined rectangular waveforms, figure 6, were used as test gust inputs. These test 
gusts were chosen because they are easily generated on the hybrid computer, have a 
wide frequency spectrum, and are repeatable. Also, transport lags between sensor sta- 
tions can be readily simulated for this form of gust. 

The general form of the feedback loops in this series was 
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The gains for die various control configurations are listed in appendix C. 

Figure 7 illustrates the effects of utilizing the different control surfaces upon the 
magnitudes of the sensed variables. The root mean square of each variable is plotted for 
each combination of active controls used on the feedback system. The values were 
obtained for the test gust disturbances with magnitudes of 1.77 m/sec. The points are 
plotted against the number of active controls with interconnecting lines added to aid 
visualization of the results. 

Since the phugoid period of 16 seconds for the basic airplane is within the 20 second 
simulation interval, the results include contributions from both long and short period 
motions. However, since a finite time interval was used in computing the performance 
information, very long period motions of the closed loop systems, which might appear, 
would not grow fast enough to be detected. An adjustment of some of the feedback gains 
might be used to stabilize these modes without increasing the gust response. 

Care should be exercised in interpreting the results presented in figure 7. The 
aircraft with controls locked at the trim positions (free aircraft case) exhibits relatively 
large magnitudes for all the variables presented. These magnitudes are due to the 
poorly damped phugoid mode being excited by the gusts. Most of the control configura- 
tions investigated show a sizable improvement relative to this case. However, a piloted 
aircraft would certainly not achieve the magnitudes presented for the free aircraft case. 

The normal acceleration magnitudes may also be misleading. In the simulation 
model, the normal acceleration has a term that is proportional to the gust velocity. 

Step changes in the gusts produce a step change in the normal acceleration. Since the 
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control actuators have lag, they cannot react immediately to these transitions. In the 
simulation several of the normal acceleration measurements are taken immediately after 
the gust transitions. Therefore, due u> the form of the gusts and the method of measure- 
ment, the measured RMS magnitude of the normal acceleration (for all control configura- 
tions) is higher than the magnitude to be expected for more realistic gust forms. For 
these reasons, the results presented in figure 7 are complemented by a closed loop pole 
analysis and a frequency response analysis in later sections of this report. 

The spoiler used alone is not as effective as flap and elevator alone cases. The 
direct lift modulation produced by the spoiler improves the normal acceleration and 
flight path angle variations with no substantial change in the other measured variables. 

The flap only control configuration shows improvement of all measured variables. 
The magnitudes of the accelerations a n and a^ are smaller than for the spoiler alone 
and elevator alone cases. Although the flap is not intended for pitch control, a substan- 
tial improvement in the pitching response is apparent. 

The elevator only configuration shows improvement of all the variables presented 
over the controls locked case and exhibits the least pitch angle and rate variation of all 
the control configurations investigated. This result is expected since the elevator is the 
obvious pitch control device; however, the lag involved in producing lift and drag changes 
is evidenced by greater accelerations than obtained by the flap only control configuration. 
The airspeed and flight path angle variations for the elevator case are similar to the 
variations for the flap case. 

Flight path angle, speed variatioi s, and normal acceleration are each improved by 
the addition of a second control device. However, these improvements are often accom- 
panied by an increase in the pitch motion. 

The control configuration with all three control devices active has less flight pitch 
angle and speed variations than any other configuration. The normal and axial accelera- 
tions and the pitch angle variations are not improved over the values obtained for two 
controls active. The pitch rate is increased over the cases with either flap or spoiler 
locked. 


Pole Locations for Various Control Configurations 

The roots of the characteristic equation (eigenvalues or poles) of the simulation 
model were examined. Figure 8 illustrates the pole locations in the complex plane for 
the various controls involving different combinations of control devices. In each case, 
any mode of the aircraft and control system which is not excited by control activity is 
deleted from the analysis. The leleted modes include the u t mode that is excited only 
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by gusts and those actuator modes for control devices which are locked in place. The 
mode shapes (eigenvectors) were also examined and the results are discussed. 

The pole locations for the aircraft with locked controls is given in figure 8(a). 

Five poles are shown; one of these is the result of the first order lag representation for 
the wing to tail transport lag of the downwash. Two poles form a very lightly damped 
oscillatory mode with period of 16.5 seconds. This corresponds to a phugoid motion but 
has a large amount of angle of attack variation a ^ ©). One pole corresponds to an 
exponentially decaying mode with time constant of 2 seconds. This is principally a 
plunging (angle of attack variation) mode. The two remaining poles form a well damped 
oscillatory mode with undamped natural frequency of 2.0 rad/sec and damping ratio 0.96. 
This is the short period mode with little speed and flight path angle variation. 

The control actuator dynamics are also represented on this figure. The elevator 
and spoiler actuators are first order lags with time constants of 0.2 second each. The 
flap actuator is a first order lag with 0.5 second time constant. 

Figures 8(b) to 8(h) illustrate the closed loop pole locations for the various control 
configurations. These locations show some similar features among the different cases. 
With the exception of the spoiler only control case, the poles form similar patterns. One 
pair, which may be either real or complex, appears near the origin. These poles may 
represent modes with time constants in the order of 100 seconds. These modes are 
primarily involved with speed and flight path angle variations as is the conventional phu- 
goid motion. Each case involving the elevator as a controlling device has a relatively 
large angle of attack variation for these modes whereas the cases of flap only control and 
flap a' u spoiler control have relatively large pitch angle variation. This shift between 
angle of attack and pitch angle variations indicates the elevator's effectiveness in reducing 
pitch variations. 

A second pair of poles appears with undamped natural frequency between 0.7 and 
1.0 radian per second and damping ratio between 0.45 and 0.82 for the cases considered. 
This oscillatory mode involves angle of attack and pitch angle variations of approximately 
equal magnitudes (and phase). The magnitude of the speed variation is approximately 
0.15 m/sec per degree of pitch angle variation. This coupling between pitch angle and 
.- aed is a bit larger than would be expected for a short period mode. 

Either one or two, real or complex poles appear in the vicinity of the -2 point on 
the real axis. These modes are characterized by essentially equal pitch angle and angle 
of attack variations, large pitch rate magnitude, and small speed and flight path angle 
variations. Flight path angle and speed variations are 0.1° and 0.6 m/sec, respectively, 
per pitch rate magnitude of one degree per second. These characteristics are typical 
for conventional short period modes. 
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The remaining modes are easily identified with actuator dynamics with little 
involvement of aircraft flight variables. 

The spoiler control case, figure 8(c), is the exception to the above general remarks. 
Two pairs of complex poles result with undamped natural frequency of approximately 
0.46 rad/sec each and damping ratios of 0.9 and 0.3. Both of these modes have relative 
speed variations of 0.36 m/sec per 1° of pitch variation. The closeness of the natural 
frequencies of the complex modes indicates that this closed loop system has very poor 
control characteristics and it would be unsuitable for use. 

Frequency Response Characteristics 

The response of the various closed loop systems to gusts were examined over a 
broad frequency spectrum. Figures 9(a) to 9(h) are frequency response plots of the 
amplitudes of normal acceleration of the center of gravity, pitch angle, flight path angle, 
and speed variation with respect to unit magnitude horizontal and vertical sinusoidal 
gusts. The plots are given for the cases of locked controls, elevator only control, ele- 
vator with flap control, elevator with spoiler control, and elevator with flap and spoiler. 
Comparison of these plots with different control configurations gives some insight into 
the effects and benefits of the utilization of the individual control devices. 

The controls locked case exhibits a high phugoid resonance peak at 0.45 rad/sec 
due to poor damping in this mode. However, there is no distinguishable short period 
mode. 

It is evident from the figures that the aircraft with locked controls is much more 
sensitive to horizontal gusts than to vertical gusts. The large horizontal gust responses 
are related to the high lift coefficient, C L , and the positive aerodynamic derivative 
C m . A head-on gust produces a large lift force due to the size of C^. Also, the 

head-on gust produces a nose-up moment. The resulting pitch-up increases angle of 

attack, increasing lift and drag supplementing the values directly attributable to the gust. 

On the other hand, a vertical, up-gust produces a nose-down moment, since C is 

m a( 

negative. Ths resulting pitch down motion decreases angle of attack and therefore lift 
and drag are diminished. 

It is also known that a positive C m has a pronounced destabilizing effect upon 

the phugoid mode. This effect is evidenced by the high phugoid peak respon *es due to 
both horizontal and vertical gusts. 

The use of elevator feedback eliminates the strong phugoid response peak. Dis- 
tinct phugoid and short period peaks appear at approximately 0.15 rad/sec and 
1.0 rad/sec, respectively. These peaks are relatively broad; indicating good damping. 
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Addition of flap or spoiler to the control configuration generally decreases the aircraft 
response to gusts. However, the pitch angle and speed variations in response to vertical 
gusts show increases in amplitude with these additions. These increases reflect the 
tradeoff activity present in the design process. 

Figures 9(a) and 9(b) illustrate the pitch angle variation due to horizontal and verti- 
cal gusts, respectively. The elevator only control law reduces the peak amplitudes 
considerable. The pitch response of this closed loop system is most sensitive to hori- 
zontal gusts at phugoid frequencies. Addition of flap or spoiler results in diminished 
phugoid response to either gust input while the short period response to horizontal gust 
is virtually unchanged. However, the short period response to vertical gusts is more 
than doubled with the addition of flap or spoiler. This increase is the result of tradeoffs 
occurring between the angular and the acceleration terms in the performance index. 
Examination of the pitch response characteristics shows no preference between the flap 
and the spoiler as active control devices. 

Figures 9(c) and 9(d) illustrate the normal acceleration response to horizontal and 
vertical gusts, respectively. Since the gusts were assumed to act upon a single, 
unswept lifting- line located at the center of gravity, i.e., at a single point in the longi- 
tudinal plane, there is no cancellation of gust due to distribution of the sinusoidal forcing 
function along the wing chord and length of the fuselage. Also, the effects of varying 
downwash angle and stream velocity in the vicinity of the tail diminish at high frequency 
due to the method of modeling the wing-to-tail transport lag. Therefore, the computed 
transfer functions for normal acceleration approach a constant level at high frequency. 

In reality, however, the normal acceleration responses trailoff towards zero. The 
mentioned deficiencies in the simulation model cause the magnitude of the normal accel- 
eration to be overestimated at high frequencies. 

The elevator only controller eliminates the large resonant peaks in the normal 
acceleration transfer functions but increases the response in the short period region. 

The addition of the flap to the control configuration reduces the responses over all fre- 
quencies. However, the spoiler has much more effect at higher frequencies. These 
results suggest the desirability of utilizing the spoilers in the control system. 

Figures 9(e) and 9(f) illustrate the flight path angle response to horizontal and ver- 
tical gusts, respectively. The elevator only control reduces the phugoid resonance peak 
to a smaller broad peak. The addition of the flap or spoiler to the control configuration 
generally reduces the flight path angle response throughout the frequency range shown, 
especially at the phugoid frequency. The spoiler has slightly more effect than the flap. 

Figures 9(g) and 9(h) illustrate the speed response transfer functions for horizontal 
and vertical gusts, respectively. The elevator only feedback law results in a broad peak 
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at the phugoid frequency and greatly reduced response to horizontal gusts at the short 
period frequency. The addition of the spoiler to the control configuration improves the 
low frequency responses by a small amount with some increase in the short period 
region. Including the flap in the control configuration results in a large reduction in 
response to horizontal gusts, especially at low frequencies. However, this reduction is 
accompanied by a substantial increase in response to vertical gusts at very low frequen- 
cies. This increase is the result of small angle approximations used in developing the 
linearized model used in this study. The effectiveness of the flap in improving speed 
response as evidenced by figure 9(g) indicates the desirability of its incorporation into the 
feedback control system. 

In summary, the frequency response to gust plots have demonstrated the effects of 
feedback control laws in reducing STOL motions in gusty conditions during the landing 
approach. The elevator provides stability as evidenced by the obvious reduction in reso- 
nant peaks. The addition of the spoiler to the control configuration is effective in 
reducing normal acceleration and flight path angle responses. The utilization o f the flap 
effectively reduces speed variations. This analysis indicates that it would be desirable 
to include elevator, flap, and spoiler in a feedback control scheme for the STOL aircraft. 

Time Histories 

Figure 10 gives the time histories of the test gust signals and the corresponding 
control surface deflections for the case of flap, spoiler, and elevator control. Close 
examination of the first few seconds of simulation reveals the nature of the control move- 
ments in response to the gusts. Vertical gusts are counteracted mainly by a spoiler 
deflection accompanied by a small retraction of the flap. These motions tend to reduce 
the wing lift coefficient in regions of updraft. Horizontal gusts are counteracted by 
simultaneous deflection of spoiler, retraction of the flap, and downward movement of the 
elevator. In regions of head-on gust, the spoiler and flap motions reduce the lift and 
drag coefficients of the wing. The elevator motion counteracts the tail £ wnload produced 
by the increase in dynamic pressure at the stabilizer. 

Figure 11 compares the time histories of several variables for the cases of locked 
controls and flap, spoiler, and elevator feedback controller. The test gust signals are 
used in these simulations. Excitation of the undamped long period mode of the open-loop 
system is readily apparent in the figure. The transitions evident in the normal accelera- 
tion result from the step changes in the horizontal and vertical gust velocities. Hie feed- 
back control law effectively reduces the pitch angle, flight path angle, and speed variation 
of the aircraft. The normal acceleration retains its sharp transitions due to the 
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rectangular nature of the gusts and the limited bandwidth of the actuators. There is, 
however, a definite reduction in acceleration magnitude immediately after these 
transitions. 


Effects of Sensor Configuration 

A series of feedback gain matrices were computed for the control configuration 
utilizing elevator, flap, and spoiler surfaces and various combinations of sensors. The 
system obtained for pitch angle and rate feedback was used as a basis for comparison. 
Normal acceleration, axial acceleration, relative airspeed at the center of gravity, and 
relative airspeed and angle of attack measured at the nose of the airframe were each 
added to the pitch angle and rate feedback system and the feedback matrices were com- 
puted. In addition, configurations utilizing both airspeed at the nose and normal accel- 
eration at the C.G. and both airspeed at the nose and angle of attack at the nose were 
used. Appendix D lists the gains for various sensor configurations. 

The placement of sensors at the nose of the aircraft necessitated changes in the 
simulation. Rather than attempt prediction of the nose location sensor signals, the time 
base for the gust signals was moved from the C.G. to the nose. Then a transport lag was 
simulated to represent the delay of the gust between the nose and the wing. These changes 
resulted in a shift in the relative phasing of the performance variable measurement times 
and the times of the transitions in the wing forces produced by the gusts. The result was 
solution control laws which had very large spoiler deflections relative to the other con- 
trols. In order to reduce the spoiler deflection magnitude to a level approximately equal 
to that of the previously discussed systems, the weight on the spoiler command signal 
squared was increased from 20 to 100. 

Figure 12 illustrates the effects of utilization of various sensor combinations 
upon the observed RMS values of several flight variables in response to the test gust 
signals. The cautions previously given relative to interpreting figure 7 apply also to 
figure 12. The addition of an airspeed sensor to the basic system results in the best 
improvement in the response of the pitch angle and inertial speed variations, and a nota- 
ble improvement in the normal acceleration and flight path angle response. The speed 
response is better when the airspeed sensor is located at the nose of the airframe, Up, 
rather than at the center of gravity, u^. These improvements are the results of a large 
airspeed to elevator gain which decreases the download on the tail induced by head-on 
gusts. This loop compensates for the large C of the basic aircraft. 

U H 

The addition of normal acceleration, a R , to the basic feedback system yields a 
good improvement in the normal acceleration response, as might be expected, with 



small improvements in pitch angle and flight path angle variations. However, with this 
sensor combination, the speed response is degraded. 

The additiou of axial acceleration, a^, feedback to the basic system results in 
small improvements in the pitch angle and speed responses. 

The angle of attack sensor placed at the nose of the airframe, dtp, contributes 
only a small improvement in the speed response relative to the basic system. Examina- 
tion of the gains for this case in table in shows that a sensed upgust results in an exten- 
sion of both the flaps and spoiler and an upward motion of the elevator. The net result 
of these control motions is an increase in lift and drag and a nose-up moment. These 
control forces are opposite to those required to cancel the effects of an up-gust. Ibis 
apparent contradiction in design specification and results is caused by the contamination 
of the vane signal by the pitch rate of the airframe. A pitch-up motion, which is pro- 
duced by a horizontal gust, results in a negative angle of attack indication by the angle of 
attack sensor. The control motions reduce the lift and drag and counter the pitch-up 
tendency of the airframe. The optimization process has resulted in a feedback system 
which uses the intended angle of attack sensor as a pitch rate sensor in this case. 

Two configurations utilizing four sensed variables were studied. The first had 
normal acceleration of the center of gravity and airspeed measured at the nose of the 
airframe in addition to the pitch angle and rate of the basic system; the second substi- 
tuted the angle of attack at the nose for the normal acceleration. The observed RMS 
response to the test gust signals for the normal acceleration and flight path angle varia- 
tions were improved for these two sensor configurations. However, the speed and pitch 
angle variations were degraded relative to the previous cases utilizing airspeed sensors. 
The former of these two cases (utilizing normal acceleration) had the best overall 
responses in normal acceleration and flight path angle. The latter (utilizing angle of 
attack) had slightly better pitch angle and speed responses. Examination of the feed- 
back gains for this last case shows that the angle of attack sensor is indeed behaving as 
an angle of attack sensor. That is, the inclusion of airspeed sensing has effectively 
removed the corrupting pitch rate signal. 

Pole Locations for Various Sensor Configurations 

Figures 13(a) to 13(d) illustrate the closed loop pole locations for each of the 
sensor configurations. The baseline case utilizing only pitch angle and rate feedback 
has an overdamped phugoid mode with 3.64 seconds time to half amplitude and short 
period damping ratio of 0.55. The addition of sensors generally decreases the undamped 
natural frequency of the short period but there is little effect on the damping ratio. The 
addition of either normal acceleration, axial acceleration, or angle of attack sensors at 
the nose of the airframe generally doubles the phugoid time to half amplitude. Each 



sensor configuration utilizing an airspeed sensor has a pole practically on the origin. 

This condition of marginal stability arises due to conflicting requirements of aircraft 
insensitivity to airspeed changes caused by horizontal gusts and speed (phugoid motion) 
stability produced by drag coefficient variation with speed. 

Frequency Responses 

Figure 14 gives frequency response curves of the magnitudes of pitch angle, nor- 
mal acceleration, glide path angle, and inertial speed variations for the various sensor 
configurations. The responses are computed using sinusoidal horizontal and vertical 
gusts of unit amplitude. These curves generally agree with the responses observed for 
the test gust disturbances and they include much lower frequencies than could be con- 
tained in the 20 second simulation interval. 

The addition of normal acceleration to the basic sensor configuration of pitch angle 
and rate resulted in a significant decrease in the normal acceleration response to both 
horizontal and vertical gusts. However, this alleviation is accompanied by increased 
responses of the pitch angle, speed, and flight path angle at the phugoid frequency. 

The addition of the angle of attack at the nose of the aircraft to the basic sensor 
configuration resulted in small improvements in the aircraft responses to horizontal 
gusts and, in general, degradation of the responses to vertical gusts. The gains in this 
case were such that a sensed upward gust causes trailing edge upward motion of the ele- 
vator. In the event of a head-on gust, the airframe pitches up, inducing a down-gust 
indication from the angle of attack sensor. The resulting dovmward motion of the elevator 
tends to cancel the pitching moment due to the gust. Therefore, this control law has 
improved response to horizontal gusts. However, updrafts cause an up- gust indication. 
The resulting upward motion of the elevator tends to pitch the nose up, This increases 
the aircraft's response to vertical gusts. 

The addition of axial acceleration to the basic sensed variables resulted primarily 
in a decrease of the pitching response to horizontal gusts, Increase of the pitching 
response to vertical gusts, and decreases of speed variations in response to both types 
of gusts. 

The sensor configuration utilizing airspeed at the nose in addition to the basic pitch 
angle and rate sensors results in improved flight path angle responses, especially in 
response to horizontal gusts. There is no apparent difference in moving this sensor to 
the center of gravity. The addition of angle of attack or normal acceleration sensors 
further improves the flight path angle responses, with the accelerometer being the more 
effective. However, very low frequency responses deteriorate with these last additions, 
indicating the approach of neutral stability. 
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The sensor configurations utilizing airspeed display considerable improvement in 
speed response to horizontal gust and varying degrees of increase in speed response to 
vertical gust. The smallest response to either gust type an obtained with the airspeed 
sensor located at the nosp. The addition of the angle of attack sensor or normal accel- 
eration sensor to this configuration results in laiger speed responses to either horizontal 
or vertical gusts. The system acquires large responses to very low frequency vertical 
gusts with these additions. 

The sensor configurations utilizing an airspeed sensor show significant reductions 
in the pitch angle response to horizontal gusts, accompanied by increases in response to 
vertical gusts. Hie least pitch variations of the sensor configurations studied occurred 
with the airspeed sensor located at the nose in addition to the pitch angle and rate feed- 
back sensors. Moving the airspeed sensor to the center of gravity nearly doubles the 
height of the response peaks at 1 rad/ sec. Hie pitch angle response is progressively 
greater for the configurations using the angle of attack at the nose, and normal accelera- 
tion in addition to the airspeed at the nose. 

The configurations utilizing an airspeed sensor are characterized by improvement 
of the normal acceleration response to horizontal gusts. Of these configurations, the 
responses are greatest for the case of airspeed sensor at the center of gravity, and are 
progressively less for the cases of airspeed sensor at the nose, airspeed and angle of 
attack at the nose, and airspeed at the nose and normal acceleration sensor. 

CONCLUSIONS 

1. The aircraft examined in this report which flies at a high lift coefficient and has 
a large pitching moment change due to changes in airspeed, was found to have large 
responses to horizontal gusts. 

2. Modeling the transport lag of the downwash variation between the wing and the 
tail with a first order lag makes the aircraft dynamics a fifth order system. For the 
aircraft studied in this report, modes similar to the conventional phugoid and short 
period modes can be identified. In addition, there is a first order plunging (angle of 
attack and speed) mode in the unaugmented aircraft dynamics. 

3. Feedback control using the elevator results in decreased aircraft response to 
gusts. This is largely due to increased stability of the phugoid mode. 

4. Inclusion of spoilers in feedback control systems decreases the normal accel- 
eration and flight path angle responses to gusts. 

5. The use of the trailing edge element of the flaps in feedback control systems 
reduces the low frequency pitch angle and speed responses to horizontal gusts. 
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6. The use of a fast response airspeed sensor in the control system reduces air* 
craft responses to horizontal gusts. 

7. The effective use of an angle of attack sensor located at the nose of the aircraft 
requires either compensation or a sufficient number of other control loops to minimize 
corruption of the measurement by the pitching motion of the aircraft. 

9. The closed loop systems obtained in this study tend to have marginal speed 
stability. This characteristic might be corrected by minor change.? in the control sys- 
tem gains. 

9. The use of test gust signals of rectangular nature in the design process has 
resulted in closed loop systems having reasonable responses to s' • •■ wdal gusts over a 
broad frequency range. 

10. The inclusion of penalties on the feedback gains in the per .mce index of an 
optimization routine prevents large gains and excessively high frequency, lightly damped 
closed loop dynamics. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., May 20, 1974. 



APPENDIX A 


AIRCRAFT SIMULATION MODEL 

A mathematical model of the externally blown flap STOL aircraft was developed 
along conventional lines with some variance in the treatment of the transport lag between 
the wing and the tail. The longitudinal dynamics of the aircraft relative to the wind axes 
coordinate system are represented by a system of differential equations: 

mVr = L - mg cos F (lift equation) 

mV = -D - mg sin I' (drag equation) 

2 2 ** 

mkyC 0 = M (pitching equation) 

The flight path angle, r, may be replaced by the equivalent expression, 0 - Aj. 

In addition, each variable may be partitioned into two parts: a nominal variable cor- 
responding to a normal flight trajectory and a relative variable being the variation of the 
original variable about its nominal value, for example, 0 » 0 Q + A0 ■ 0 q + 0. The 
aerodynamic forces may also be partitioned into portions contributed by the wing and by 
the tail, for example L = L w + L t> The incorporation of these items in equation (1) 
yields an alternate form for the aircraft dynamical equations. 

m(V 0 + v)(f 0 + 0 - dj) - + L t - mg cos(r o +0 - aj) 

m(v 0 + v) a -D w - D t - mg sin(r o + 0 - 

mk2e 2 (0 o + 0) - M w + M t (A2) 

The lift, drag, and moment produced by the wing ^L w ,D w ,M w j are functions of 
the angle of attack of the wing, the relative stream velocity at the wing, and the control 
surface deflections. The angle of attack is composed of the inertial angle of attack, 
caused by the motion of the aircraft through space, and contributions induced by gusts. 
The relative stream velocity at the ving is composed of the inertial velocity, V, and 
contributions induced by gusts. The lift produced by the wings may be expressed, in 
general, as 

L w * L w 0 * i w (V.V v .V H )SC Lw (V,a I ,V v ,V H ,e.,M) (A3) 
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Expansion of this expression by Taylor series and discarding terms above the first order 
results in the following linear expression for the lift of the wing: 


L = L + i pv£ S 
w w Q 2 ° 


(^Ly, + C L ,U + C L «I + Ct Uy + Ct U H 

\ ^0 V «W U V U H 

w n w 


where 


+ Ct 5S + Ct 6f 

l 68v, L a w 


V AV V v 

*7' V" 11 ' *o ’ 


(A4) 


Similar linearized expressions can be derived for the drag and moment of the wing. 

The lift, drag, and moment produced by the tail are functions of the 

relative stream velocity and direction at the tail and the elevator deflection. The angle of 
attack of the horizontal tail is composed of the angle of attack of the wing, the downwash 
angle including gust effects, and the pitch rate of the aircraft. The relative stream 
velocity at the tail is composed of the inertial velocity of the aircraft and gust contribu- 
tions. The effect of downwash variations reach the tail at a time (r = f t c/V 0 ) after the 
variations are caused at the wing. Also, the stream velocity variations at the tail due to 
the gusts occur at a time, r, after the gusts effect the wing. The lilt produced by the 
util may be expressed, in general, as 

L t (t) = - r),V H (t - r)jsc Lt (a I (t),e(t),6(t),6e(t)) (A5) 

Expansion of this expression via Taylor series and discarding terms above the first 
order results in the following linear expression for the lift of the tail: 


L t = \ 


+ 2 pV 5 S !* C L t ( u u t) + % (“I - e ) + c hj 

l 0 t 0 


6 + c. 


J 6e 


Se 


(A6) 


where 


8u f 9Uf 

u t (t)a d4 u H (t - T) + a^ u v (t - T) 

The variable u t is the portion of the stream velocity at the tail due to the gusts. 
Similar linearized expressions can be derived for the drag and moment of the tel 
(D. and MA 

\ t w 
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APPENDIX A - Continued 


The downwash angle is a function of the angle of attack of the wing, the gust veloc- 
ities, the aircraft velocity, the spoiler deflection, and the flap deflection. 


€ (t) = H«!(t - T) + ft u(t - T) + ^u v (t - T) + ^-U H (t - T) 

J. -5l. fisft - rl 4- M(t . r \ 


(A8) 


The linearized expressions for the lift, drag, and moment components (eqs. (A4) 
and (A6)) are substituted into the aircraft dynamic equations (eq. (A2)). Small angle 
approximations are made for the trigonometric functions. For the steady flight condition 
the following simplifying relations are utilized: 


V 0 
V° 
© 0 = 0 






(A9) 


L o ■ \ * \ ■ cos r „ 

D o * D w 0 + \ = -“B sln r o 

M = M +M; =0 

0 0 > 

Second order terms are dropped and the resulting linearized equations may be written as 


-«i = ^-(e - a^ sin r o * © + 1 m v o s / c L a a i + c l . 6 
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APPENDIX A - Concluded 


is '2™ V ° S f 1 D a w “l + + [" C D Uw + 2C D> + C D U V + C D U H 

L *- “* V w Hw 

1 C »«.«‘ * C D** * C D 58w » s * 2C D t “t - C D at ( f * - “l>° 8 r 0 (AlOc) 

The transport lags in equations (A7) and (A8) are treated differently from the usual 
procedure (see ref. 4). The lags are approximated by first order lags whose time con- 
stants are equal to the time of delay of the transport lags. The relations used to approx- 
imate e and u t are 


= fe + -It- + 4t- «u + •§£ u + 6f + 


Cl. . 

^ U t = 


^“I + ^ U V + ^ U H + ^ 


9u t 9u t l 
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(All) 



APPENDIX B.- AIRCRAFT PARAMETERS 



= 3.203 m (10.51 ft) 

= 74.45 m 2 (801.36 ft 2 ) 

* 23.69 m 2 (255.0 ft 2 ) 

= 3.2 

= 1.6 

* 25 022 kg (1715 slugs) 

* 1.31 

= 1.225 kg/m 3 (0.002378 slug/ft 3 ) 
= 9.805 m/s 2 (32.17 ft/s 2 ) 

= 60° 

= 10 ° 

= 1.87 
= 15° 

= 13.89° 

= -25° 

= 5.16° 

= 4.7633 

* 0.3231 

* -0.5160 
= 0.4826 
= -0.0326 
= 0.5160 
= 4.2807 

* 0.2905 
= -3.88° 

* 6.12° 

= 35.41 m/s (116.17 ft/s) 
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= 5.7 rad" 1 

= 4.0 rad" 1 

= 1.7 rad" 1 

= -3.50 

= 1.00 

= 1.13 

= 2.8 rad" 1 

= 3.5 rad" 1 

= -0.6 rad" 1 

= -2.89 rad -1 

= -0.408 rad" 1 

= 1.581 rad' 1 

= 0.5590 
= -0.3097 rad 
= 0.2746 
= -0.2834 
= 1.6146 rad" 1 

* 0.5185 rad" 1 

* -5.3700 rad" 1 

= 0.4166 (rad/sec)" 1 
= 0.1014 (rad/sec)" 1 
= -1.3252 (rad/sec)" 1 
= 0.3062 rad -1 
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6e 
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6e 


X 

X 
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9u t 

dc 

de 
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= 0.0687 rad * 

= -0.9783 rad" 1 
= 5.6050 

* 1.6942 

= -0.0175 

= 6.4174 

= -0.6501 

= 1.7027 

= 0.9977 

= 0.0677 

= -0.0098 

= -0.4425 

= 0.2 sec 
= 0.2 sec 
= 0.5 sec 


APPENDIX C. - GAINS FOR VARIOUS CONTROL CONFIGURATIONS 


Elevator Only Controller 

oe c = 0.1979 a n - 0.5115 6 - 0.0163 0 + 1.2642 u A - 0.1878 a x 
Spoiler Only Controller 

6s c = 1.5012 a n + 0.0220 0 - 0.2949 0 + 0.6678 u A - 0.2343 a x 
Flap Only Controller 

6f c » -0.0208 a n - 0.1035 0 - 0.56980 0 - 0.8494 u A - 0.4984 a x 
Elevator and Spoiler Controller 

5e n = 0.2391 a - 0.0374 0 - 0.2849 0 + 1.2540 u. + 0.3396 a v 
c n AX 

6s c = 0.7237 a n - 0.1142 0 + 0.2702 0 + 0.5574 u A + 0.2592 a x 
Elevator and Flap Controller 

6e c = -0.0491 a n - 0.0738 6 - 0.1001 0 + 1.0275 u A + 0.0387 a x 

6f c = -0.1492 a n + 0.2446 0 - 0.8521 0 + 0.4756 u A - 0.2899 a x 
Spoiler and Flap Control ler 

6s c = 0.7125 a n - 0.5137 0 + 0.5018 0 - 0.0033 u A + 0.2706 a x 

6f c = 0.0553 a n - 0.2910 O - 0.4811 0 - 0.8547 u A - 0.4367 a x 
Elevator, Spoiler, and Flap Controller 

6e c = 0.1949 a n + 0.2488 0 + 0.2135 0 + 1.0216 u A - 0.2143 a x 

6s c = 0.6300 a p - 0.0989 0 + 0.2478 0 + 0.3660 u A - 0.3990 a x 

6f c = -0.0251 a n + 0.2386 6 - 0.5350 0 - 0.4027 u A - 0.0889 a x 
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APPENDIX D. - GAINS FOR VARIOUS SENSOR CONFIGURATIONS 


4 

1 
js 

i 

Pitch Angle and Rate Only 

6f = -0.8698 0 - 1.0033 0 

6s c = 0.7633 0 + 7408 0 

6e c = 1.7368 0 + 0.7324 0 

With Normal Acceleration 

6f = -0.2938 a„ - 0.8821 © - 1.0025 0 
c n 

6s„ = 0.9888 a - 0.3343 0 + 0.2485 © 
c n 

6e c = 0.8505 a h - 1.3370 0 - 0.1589 0 

With Airspeed at C.G. 

6f c = -0.0340 0 - 0.6664 0 - 0.4253 u A 

6s c = 0.5283© + 0.3308 0 + 0.3893 u A 

6e c = -0.4929 0 + 0.1900 0 + 1.0931 u A 

With Axial Acceleration 

6f„ = -0.4759 0 - 0.9072 0 + 0.4755 a v 
c X 

6s. = 0.7183 0 + 0.6706 © + 0.0258 a v 
c X 

6e„ = 1.0721 0 + 0.3726 0 - 1.8043 a v 
c X 
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APPENDIX D - Concluded 


With Airspeed at Nose 

6f c = -0.4153 0 - 0.5166 0 - 0.4282 Up 

6s„ = 0.5793 © + 0.1827 0 + 0.2908 u- 
c F 

6e„ = -0.0866 0 - 0.0250 0 + 1.0320 u- 
c r 

With Angle of Attack at Nose 

6f = -0.8277 0 - 1.0057 0 + 0.3754 a v 
c r 

6s c = 0.6282 0 + 0.5092 0 + 0.211 f a p 

6e c = 1.3419 0 + 1.0476 0 - 0.6551 a p 

With Airspeed at Nose and Normal Acceleration 

6f = -0.2575 a„ - 0.2421 0 - 0.3619 0 - 0.3550 u„ 
c n r 

6s _ = 0.8365 a - 0.0501 9 + 0.2024 0 + 0.1678 u F 
c n r 

6e c * -0.0652 a R + 0.2805 0 - 0.1670 0 + 1.0908 u p 
With Airspeed and Angle of Attack at Nose 

6f c = -0.4075 0 - 0.4012 0 - 0.5519 u p - 0.2200 <* F 

6s c = 0.6090 G + 0.1426 0 + 0.3553 u p + 0.3093 

6e c = 0.0029 0 - 0.1102 0 + 1.0436 u p + 0.0226 a p 
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Wing spoilers 



Figure 1.- Three-view drawing of simulated airplane. All linear dimensions 

are in meters (feet). 



Figure 3.- Simulation variables. 


f 


1 



Figure 4.- Conjugate gradient algorithm logic flow chart. 
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(a) Normal acceleration. (b) Axial acceleration. 

Figure 7.- Observed variable magnitudes for different control configurations. 
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and spoiler 
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and spoiler 


1 2 3 0 1 2 

No. of active controls No * of active controls 

(c) Airspeed variation. (d) Flight path angle variation. 

Figure 7.- Continued. 
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(e) Pitch rate. (f) Pitch angle variation. 

Figure 7.- Concluded. 
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Figure 8.- Pole 
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locations for various control combinations. 







Figure 8.- Continued. 







Frequency, rad/sec 

(a) Pitch angle response to horizontal gust. 

Figure 9.- Frequency responses to unit amplitude sinusoidal gusts 
for various control configurations. 
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Figure 9.- Continued. 
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Figure 10.- Test gusts and control deflections. 
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Figure 11.- Feedback system effectiveness. 
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(a) Normal acceleration. (b) Pitch angle. 

Figure 12.- Addition of sensors to pitch angle and rate feedback control system. 
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(c) Flight path angle. 


(d) Speed. 


Figure 12.- Concluded. 
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(f) Pitch angle, pitch rate, and angle of attack at nose. 


Figure 13.- Continued. 
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Frequency, rad/sec 

(a) Pitch angle response to horizontal gust. 

Figure 14,- Frequency responses to unit amplitude sinusoidal gusts 
for various sensor configurations. 
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— ■- — ■ ■ Pitch angle and rate feedback 

— - With angle of attack 

— — With airspeed at nose 

With angle of attack and 
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(a) Pitch angle response to horizoultl gust - Concluded. 
Figure 14.- Continued. 
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(c) Normal acceleration response to horizontal gust. 
Figure 14.- Continued. 
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(c) Normal acceleration response to horizontal gust - Concluded. 
Figure 14.- Continued. 
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(e) Flight path angle response to vertical gust - Concluded. 
Figure 14.- Continued. 
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(g) Speed response to horizontal gust - Concluded. 
Figure 14.- Continued. 
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Pitch angle and rate feedback 
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(h) Speed response to vertical gust. 


Figure 14.- Continued. 
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(h) Speed response to vertical gust - Concluded. 


Figure 14.- Concluded. 
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